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We report the results ol exact diagonalization studies of Hubbard models on a 4 x 4 square 
lattice with periodic boundary conditions and various degrees and patterns of inhomogeneity, which 
are represented by inequivalent hopping integrals t and t' . We focus primarily on two patterns, 
the checkerboard and the striped cases, for a large range of values of the on-site repulsion U and 
doped hole concentration, x. We present evidence that superconductivity is strongest for U of 
order the bandwidth, and intermediate inhomogeneity, < t' < t. The maximum value of the 
"pair-binding energy" we have found with purely repulsive interactions is A p f, = 0.32t for the 
checkerboard Hubbard model with U — 8t and t' — 0.5f. Moreover, for near optimal values, our 
results are insensitive to changes in boundary conditions, suggesting that the correlation length is 
sufficiently short that finite size effects are already unimportant. 
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The relatively large energy scales and short coherence 
lengths involved in high temperature superconductivity 
(HTC) imply that theories of the "mechanism" must in- 
volve different considerations than the conventional BCS 
theory of low temperature superconductivity (LTC). The 
theory of LTC can be treated in the context of Fermi liq- 
uid theory, in which the strong effects of electron-electron 
repulsions are resolved at high energy, so that pairing is a 
Fermi surface instability triggered by weak, retarded, in- 
duced attractive interactions. The theory of HTC must 
treat the strong local repulsions between electrons di- 
rectly, as Fermi liquid theory is certainly not valid at 
short distances and high energies. Conversely, since the 
physics of HTC is relatively local, numerical studies of 
finite size systems can, plausibly, resolve questions con- 
cerning the mechanism so long as the system sizes are 
large compared to the coherence length, a condition that 
could not remotely be envisaged for LTC. 

In the present paper, we report exact diagonalization 
studies of the low energy states of the Hubbard model on 
a square lattice (Fig. [I}, with hopping matrix elements, 
tij, between pairs of sites, i and j, and with a short-range 
repulsive interaction, Uj, between two electrons on the 
same site. We propose to address the following sharply 
posed question, related to the physics of the mechanism: 
What form of the Hamiltonian (i.e. values of {tij} and 
{Ui}) maximizes T c subject to the constraints that 1) tij 
is short-ranged and bounded, i.e. |£y| < t for all ij, and 
that 2) the interactions are repulsive, U > for all il 
Constraint 2 represents a theoretical prejudice that HTC 
derives directly from the strong repulsive interactions be- 
tween electrons, and could be relaxed in future studies. 
The requirement that t^ is bounded avoids the trivial 
answer that, for any Hamiltonian that supports super- 
conductivity, T c can be doubled by simply doubling the 
Hamiltonian. 

The largest systems we can readily study are 4x4. 
For such small systems, there is no direct way to ex- 
tract a T c . We have thus introduced other benchmarks 
of the strength of the superconductivity that can be read- 
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FIG. 1: (Color online) Schematic representation of the inhomoge- 
neous Hubbard model: (a) the checkerboard version and (b) the 
striped version. The hopping matrix elements are t on the solid 
"strong" bonds and t' < t on the dotted "weak" bonds. The three 
labeled bonds, ab, cd, and ef are the focus of our study of pairing 
correlations. 



ily computed on finite size systems, especially the "pair 
binding energy," defined in Eq. ([3]), and the strength of 
the pair-field correlations defined in Eqs. (H))-©. 

The full optimization problem we have proposed would 
be prohibitively time consuming. We have therefore con- 
centrated on a restricted, highly symmetric subset of 
all possible Hamiltonians, with particular focus on the 
symmetric checkerboard and stripe patterns shown in 
Figs, nja) and Hps). We take the hopping matrix ele- 
ments to be t on the solid ("strong") bonds in the figure, 
and t' < t on the dotted ("weak") bonds. For t' = t, 
the system is the homogeneous Hubbard model, while 
for t' = 0, the system consists of disconnected Hubbard 
squares or ladders. Thus, as we vary t' from to t, we 
vary the "degree of inhomogeneity" . 

For instance, as shown in the contour plot in Fig. [21 the 
pair-binding energy of the checkerboard lattice with pe- 
riodic boundary conditions (PBC) is largest for U = 8t 
and t' = 0.5t. The concentration, x, of "doped holes" 
per site (i.e. the deviation from one electron per site) 
can only take discrete values; among those, the optimal 



pair-binding is largest when x — 1/16 (as in the con- 
tour plot in Fig. Eta)), somewhat smaller when x = 3/16 
(Fig. [2(b)), and is always negative (pair- repulsion) for 
x = 5/16. Indeed, of all forms of the Hamiltonian we 
have explored to date, the checkerboard Hubbard model 
with these parameters has the largest pair-binding energy 
we have found. Moreover, by changing boundary condi- 
tions to twisted PBC, shown as the solid (red) triangles 
in Fig. [SKa), we see that at fixed U = 8t, as a function of 
t' in the range 0.8 >t'/t> 0, the results are sensitive to 
change of boundary conditions only at the 20% level. 

We therefore infer that, the checkerboard Hubbard 
model in the thermodynamic limit, has a maximum value 
of the pair-binding A p b f=a i/3 for U ks 8t, t' sa t/2, and 
x ks 1/16. From analysis of the ground-state symme- 
try and of the pair-field correlations, we can identify the 
dominant superconducting correlations on this system as 
d-wave (d x 2_ y 2). 

I. THE HAMILTONIAN 

The inhomogeneous Hubbard model we have studied 
is described by the Hamiltonian 

where indicates nearest-neighbor sites, c\ a creates 
an electron on site i with spin polarization a = ±1, and 
= cl a Ci. a . The usual (homogeneous) limit of this 
model is obtained by taking £y = t and Ui = U . 

Although we have studied a wider variety of patterns, 
to be concrete we will primarily focus on two inhomoge- 
neous patterns of the hopping amplitudes: the checker- 
board lattice and the striped lattice shown in Fig. [1] Un- 
less otherwise stated, our discussion will focus on the 
case in which PBCs are applied in both the x and y di- 
rections. However, finite size effects will be estimated by 
comparing these results to those with "twisted" periodic 
boundary conditions. Specifically, for a 4 x 4 system, 
PBC means identifying the sites (n + 4, to) = (n, m) 
and (n, to + 4) = (n,m), while twisted PBC in the 
y direction means identifying (n + 4, to) = (n, to) and 
(n, to + 4) = (n + 2, to). We have also obtained results 
with open boundary conditions, but because of the large 
surface to volume ratio of the small systems studied, the 
proper interpretation of these results is unclear, and so 
we do not report them, here. 

II. RESULTS 

The exact diagonalization is performed using the Lanc- 
zos method^, which has been employed by other people 
to extract some ground state properties of a 4 x 4 (ho- 
mogeneous) Hubbard mode£^. By taking advantage of 
conservation of the z component of the spin and the C^ v 
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(a)A pf) (l/l6), checkerboard lattice 
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(b) Apj, (3/16), checkerboard lattice 
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(c)A p f,(l/16), striped lattice 



FIG. 2: (Color online) Contour plots of the pair-binding energy 
as a function of U and t' on two types of lattices with periodic 
boundary conditions. Note change of scale in c. 
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symmetry of the checkerboard and D2h symmetry of the 
striped model, we successfully reduce the dimension of 
the Hilbcrt space from around 10 8 to 10 7 states and the 
ground state energies we found for t'/t — ► 1 agree with 
those of G. Fano et al£. A complete group theoretical 
analysis which facilitates efficient implementation of the 
Lanczos algorithm can be found in Ref . 0- 

We present representative results on the pair-binding 
energy and the pair-field correlations. We have much 
more extensive tables of results for various values of U/t, 
t'/t and x, and for various choices of boundary condi- 
tions. These additional results are available from the 
authors upon request. 

A. The Pair-Binding Energy 

To better understand the pairing phenomena arising 
from repulsive interactions, we define the pair-binding 
energy: 

A pb (x) = 2E (M) - [E {M + 1) + E (M - 1)], (2) 

where N — 16 is the number of sites in the system, 
Eq(M) is the ground state energy with N — M electrons 
(i.e. M holes doped into a "neutral" half-filled lattice), 
and x = M/N is the "concentration of doped holes." We 
will focus on the case in which M is odd. Thus, a positive 
pair-binding energy means that, given two isolated clus- 
ters with a mean doped hole density x, it is energetically 
preferable to "pair" the doped holes so that one cluster 
has M + 1 and the other cluster has M — 1 doped holes. 

For a superconducting system, A p (, — > 2A m i n in the 
limit N — ► oo, where A m i n is the minimum value of the 
superconducting gap. While gaps can occur for other 
reasons (e.g. CDW formation), a superconducting state, 
as far as we know, is the unique state that generically 
produces a gap for a non-zero range of x in more than 
ID. In a superconducting state with gapless (nodal) ex- 
citations, Apt vanishes as N — > oo, but only relatively 
slowly, in proportion to AqN^ 1 / 2 in 2D, so A p (, com- 
puted on finite size systems remains a good diagnostic 
of superconductivity. Under generic circumstances in 
non-superconducting systems, the repulsive interactions 
between quasiparticles implies that A p b is negative and 
A pb ~ -TV" 1 as N -> oo. 

A positive pair-binding energy on a small system could 
also indicate a tendency for phase separation. Unambigu- 
ously distinguishing gap formation from phase separation 
can only be done by appropriate finite size scaling^, which 
is beyond the reach of the present calculations. However, 
a gross test for phase separation is possible by testing 
whether further agglomeration of doped holes is favored. 
Specifically, for M = Nx even, we define 

k n = [E (M + 2) + E (M - 2) - 2E (M)]/2. (3) 

This is a crude approximation of the compressibility^, 
which is negative in a system with a sufficiently strong 
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FIG. 3: (Color online) The pair-binding energy, A p f,(l/16), as 
a function of t' at U = 8t on a (a) checkerboard lattice and a 
(b) striped lattice with various boundary conditions. Squares rep- 
resent data with PBC in both directions; triangles represent data 
with PBC in the x direction and "twisted" PBC in the y direc- 
tion. [Inset: A closer look at A p( ,(l/16) as t' goes to zero on 
the checkerboard lattice. Notice that A p j,(l/16) becomes negative 
when t' < 0.025i.] 

tendency to phase separation. For the interesting range 
of t' and U, we never find a negative value of kn and 
hence our system is unlikely to be phase separated. 

We have computed the pair-binding energy A p f, as a 
function of t', U, and x on both the checkerboard and 
striped lattices with PBCs and twisted PBCs. Con- 
tour plots of A p b(x = 1/16) and A p b(x = 3/16) for 
the checkerboard lattice are shown in respectively, in 
Figs. (a) and (b), respectively. The global maximum 
can clearly be seen for U — 8t and t' = Q.5t. 

Apb(x = 1/16) is shown at fixed U = 8t as a func- 
tion of t' in Fig. [3] with both PBC (squares) and twisted 
PBCs (triangles) . A remarkable degree of insensitivity to 
boundary conditions is apparent for t'/t < 0.8. On the 
coarse scale of t' in the main figure, it appears that A p & 
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FIG. 4: (Color online) The doping dependence of the pair-binding 
energy, A p j,(x), as a function of t' on the checkerboard lattice. 

is positive at small t'/t for all U. However, this hides a 
subtlety at small t'/t, as shown in the inset to the figure 
in which the regime of small t' is shown on an expanded 
scale. Specifically, as shown in Fig. the pair-binding 
energy on an isolated square changes from positive to 
negative at U = U c w 4.6t. Moreover, as discussed in 
Ref. H, it follows that the pair binding is positive with 
a non-zero limit as t'/t — » for < U < U c , while 
A p b ~ — 0([t'} 2 /t) is negative and vanishes as t'/t — > 
for U > U c . 

In Fig. Apb(x) for the checkerboard model is shown 
at fixed U/t = 8 as a function of t'/t for different values 
of x = 1/16, 3/16, and 5/16. 

A contour plot of A p i,(x = 1/16) for the striped lat- 
tice is shown in in Fig. [5] (c), and A p g(x = 1/16) for 
fixed U = 8t is shown as a function of t'/t for PBC and 
twisted PBC in Fig. [3] (b). Here, the results are appar- 
ently more sensitive to boundary conditions, so inferences 
concerning the thermodynamic limit are more difficult to 
reach. However, here, too, a global maximum of A p b is 
reached for U = lit and t' = 0.7 't. In the limit of vanish- 
ing t', as follows from the results for an isolated ladder, 
Fig. O (dashed curve), A p & is positive and non- vanishing 
for < U < U c i « 3.5i, negative and order (t') 2 for 
U c i < U < U c 2 ~ 7t, again positive and non-vanishing 
for U C 2 < U < U C i ~ 15£, and finally negative and order 
(t 1 ) 2 for U > U c3 . 

B. Pair-field correlations 

We have also studied the equal-time pair-field pair-field 
correlation function defined as 

D(ij,kl) - (Aj.A w ), (4) 

where the pair-field is 

4 = 4(4 t 4i+4t4i)- (5) 
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FIG. 5: (Color online) The pair-binding energy, A p j,(l/4), as a 
function of U for an isolated 2x2 plaquette (solid curve) and a 
4x2 ladder (dashed curve). Note that in the U = oo limit one can 
show A p i,(l/4) < for both clusters^. 

ij represents the bond between a pair of nearest-neighbor 
sites, on the lattice. We focus on this correlation 
function for the largest possible separations, given the 
small system size: between a pair of parallel strong bonds 
separated by a distance 2, i.e. bonds ab and ef in Fig. [U 
and a pair of perpendicular strong bonds separated by 
a distance 3V2/2, i.e. bonds ab and cd in Fig. [T] We 
have computed these pairing correlations for both the 
checkerboard and striped models with M = 0, 2, and 4, 
i.e. for doped hole concentration x = 0, 1/8, and 1/4, 
and for a range of U/t and t'/t. 

In Fig.^a), we show the pair correlation function for 
the checkerboard lattice at fixed U/t = 8 func- 
tion of t'/t for x = 1/8 (triangles) and 1/4 (squares). 
D(ab, cd) is represented by the dashed lines in the figure, 
and D(ab, ef) by the solid lines. The same quantities are 
shown for the striped lattice in Fig. [6jb). Qualitatively 
similar results for both lattices have been obtained for 
x = 0, but the mangnitude of D is an order of magnitude 
smaller than for x > 0, consistent with the expected Mott 
insulating character of the undoped system. The positive 
sign of D(ab, ef) and the negative sign of D(ab, cd) are 
consistent with the d-wave character of the pairing; in the 
thermodynamic, at large spatial separations, the sign of 
D is determined entirely by the symmetry of the order 
parameter. 

For the checkerboard model, it is apparent that the 
magnitude of the pairing correlations are relatively weak 
both in the limit of strong inhomogeneity (t'/t <C 1) 
and of vanishing inhomogeneity (t'/t ~ 1). However, the 
strength of these correlations is only weakly dependent on 
t'/t for a broad range of intermediate values. Moreover, 
in this intermediate regime, the strength of the pairing 
is quite similar for x = 1/8 and x = 1/4. 

It is more difficult to make a clear qualitative state- 
ment about the behavior of the striped model. One 
seemingly puzzling feature of the striped results deserves 
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FIG. 6: (Color online) Pair-field pair-field correlation functions 
at U = 8t on the (a) checkerboard and (b) striped lattices. Solid 
curves represent D(ab,ef) and dashed represent D(ab,cd); trian- 
gles are data points for M = 2 doped holes, and squares are data 
points for M = 4. Note that in the striped case D(ab,cd) have 
been multiplied by ten for comparison. 



comment. In the limit t' — > 0, where the two ladders 
decouple, one might expect the correlation of the inter- 
ladder pair fields [D(ab, cd)] to vanish. However, since for 
U = 8, the pair binding energy on a 4 x 2 ladder (cube) is 
positive, as shown in Fig. [S] (dashed curve), for t' = and 
M = 2, there are two degenerate ground states, with the 
hole-pair on one or the other disconnected ladder. Thus, 
even as t' tends to 0, the ground-state wave function is a 
coherent superposition of these two states, and hence has 
substantial pair-field correlations. Were one to study the 
finite temperature properties of this system, the coher- 
ence would be lost above a relatively small temperature 
T co h ~ t'. This illustrates the dangers of uncritically ac- 
cepting evidence of strong superconducting correlations 
from the pair-field correlation function. 

We have also studied the expectation value of the pair- 
field operator between ground-states with M = Nx and 
M - 2 doped holes: 




FIG. 7: (Color online) The square of the expectation value of the 
pair annihilation operator, |A a | 2 , as a function oft' on the checker- 
board lattice. The subscript a refers to the contribution from the 
strong and weak bonds, respectively. The open symbols indicate 
the data obtained at t' = 0.99t. As in Eq. llOt . 1/e times this 
quantity is the single-mode-approximation to the superconducting 
susceptibility. 



(Ay) = (M;0|Ay|M-2;0) 



(6) 



where |M;0) is the ground state with M doped holes. 
This is not a gauge invariant quantity, in that there is 
an arbitrary choice of an overall phase; by choosing the 
wave functions real, this is reduced to an overall sign am- 
biguity. However, the internal symmetry of the pairing 
state is manifest in this quantity. 



For the checkerboard case, for all U and t' < t, the 
ground state symmetry as a function of M alternates, 
A 1 for M = 0, B x for M = 2, A 1 for M = 4. Thus, 
the pair-field operator that connects any two of these 
states must have precisely B\ {i.e. d x i_ y i) symmetry^. 
In Fig. for fixed U = 8t, we plot, as a function of 
t'/t, the gauge-invariant quantities, |A S | 2 (squares) and 
|A„,| 2 (triangles), where A s and A w are, respectively, 
the expectation value of Ay for strong bonds ij (within 
a square) and weak bonds (connecting two squares) [see 
Eq. (8)]. Again, there is a clear indication that super- 
conducting correlations are strongest for t' w t/2. 



Two peculiarities arc worth mentioning here: Firstly, 
it is puzzling that the difference in the pair amplitudes on 
the weak and strong bonds is so large, |A^|/|A S | ~ 10~ 7 , 
as clearly shown in Fig. [7] This is a qualitative point 
that warrants further study. Secondly, even in the limit 
t' — > t, \A W \ 7^ |A S |. This is a consequence of the fact 
the ground-state is three-fold degenerate^ when t' = t 
and M = 2; this degeneracy is lifted whenever t' ^ t. 
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FIG. 8: (Color online) A speculative phase diagram of the checker- 
board model at x = 1/16. 



C. Superconducting Susceptibility 

The T = superconducting susceptibility is express- 
ible as 



X 



a(M) = E 



|(M;0|A( a )|M-2;a)| 2 
E a {M - 2) - E {M) - 2n 



(7) 



where E a {M) is the energy of the ath excited state with 
M doped holes, and where we take 

1 



El/. 



(8) 



Here, the chemical potential is appropriate to the case 
in which the ground-state has M doped holes, i.e. 
[E (M - 2) - Eq(M)] > 2/i > [E (M) - E {M + 2)]. 
Because of the d-wave symmetry, we always take /y to 
be positive on vertical bonds and negative on horizontal 
bonds. We consider three possible susceptibilities: \s in 
which where, \fij\ = 1 on the strong bonds and on the 
weak, Xw in which = 1 on the weak bonds and 
on the strong, and \t in which \fij\ = 1 on all nearest- 
neighbor bonds. For simplicity, we have defined the sus- 
ceptibility with respect to adding two- holes. Another 
(more conventional) definition includes, as well, terms 
which remove two holes. 

As mentioned before, computing all the excited states 
that enter this sum is not feasible, but a lower bound esti- 
mate can be readily obtained in the "single mode approx- 
imation" , by approximating the sum by the single term 
involving the ground state with N — AI + 2 electrons. 
We refer to this as Xa MA \ In terms of the pair-field 
expectation values discussed in the previous section, 



(SMA) = |As|2/£ 

= IA,„| 2 / e , and 
A w \ 2 /2e. 



X 



(9) 



(SMA) , * 
Xt = \ A 



where e = [E (M)-E (M +2)+2fi). Thus, the quantities 
plotted in Fig. [7] can be viewed (up to a factor of 1/e) as 
the SMA to the susceptibility at constant e. Moreover, 
it is worth mentioning that they are almost exact to the 
lower bound susceptibility as long as e — ► by tuning 
chemical potential properly. 



III. INTERPRETATION 

We now discuss the implications of our results. In par- 
ticular, we are interested in using the present information 
to infer, as much as possible, the phase diagram of the 
checkerboard and striped Hubbard models in the ther- 
modynamic limit, N — > oo. The results of our analysis 
and some additional speculations lead us to propose the 
qualitative phase diagram shown in Fig. [8l 



A. The Checkerboard Model 

There are two limits in which the checkerboard model 
simplifies: U -C *', where it can be studied using conven- 
tional diagramatic methods, and t' <C t, where it reduces 
to weakly coupled squares which can be treated^ using 
degenerate perturbation theory in *' /t. This allows us to 
deduce the solid portions shown in the phase diagram in 
Fig. [5] without resorting to the present numerical results. 

For f7 = 0, and in the thermodynamic limit, the non- 
interacting Fermi-surface depends qualitatively on t'/t. 
There are four bands since each unit cell contains four 
sites and (2n/2a,0) and (0, 2it/2a) are the basis vectors 
of the reciprocal lattice. For t > *' > t' c (x), there are two 
electron pockets enclosing the M points (ir/2a, 0) and 
(0, it /2a) respectively, and one hole pocket, enclosing the 
"nodal point" (7r/2a, w/2a), as shown in Fig. [9ja). The 
hole pocket shrinks to a point as t' approaches the critical 
value for a Lifshitz transition, t' c (x). For < t' < t' c (x), 
only the two electron pockets enclosing the M points 
remain, as shown in Fig. [UJb). For x = 1/16,1/8,3/16, 
and 1/4, t' c is, respectively, 0.95t, 0.89t, 0.82*, and 0.75*. 

This Lifshitz transition appears in the conjectural 
phase diagram in Fig. [5] as the boundary between two 
Fermi liquid phases - FL1 (with two electron pockets plus 
one hole pocket) and FL2 (with only two electron pock- 
ets). In sketching this figure, we have assumed that the 
Fermi liquid phases are stable in the presence of a small 
repulsive U; it is likely that this is not strictly the case, 
since there is probably a Kohn-Luttinger instability of 
any Fermi liquid^, but if this occurs, it is on such a low 
energy and temperature scale that it can be neglected for 
present purposes. 

The small *' portion of thephase diagram in Fig. [H] was 
derived previously in Refs. U and 0. For < U <U C ~ 
4.6*, where an isolated Hubbard square has a positive 
pair-binding energy, there exists a nodeless d(2.2_ y 2)-wave 
superconducting phase. For U c < U < Ut ~ 18.6*, there 
is a third Fermi liquid phase, FL3, which has the same 
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FIG. 9: (Color online) A sketch of the Fermi surface with x = 3/16 
for (a) t' = 0.95t, showing the two electron and one hole Fermi 
pockets, and (b) t' = 0.70t, with only the two electron pockets. 



Fermi surface topology as FL2. Finally, for Ut < U 
(where the isolated square with one doped hole has a fully 
polarized spin 3/2 ground state), the system exhibits an 
exotic, spin 3/2 Fermi liquid phase, FL4, if x < 1/8, while 
it phase separates into two insulating antiferromagnetic 
charge ordered phases ifl/8<x<l/5. 

The superconducting state at small t' arises from the 
geometry of the square and the strong correlations pro- 
duced by U. However, across the phase boundary be- 
tween the superconducting phase and FL2 or FL3, it is 
reasonable to view the onset of superconductivity as a 
BCS-like Fermi surface instability. In this limit, the fact 
that the d-wave superconductor is nodcless can be un- 
derstood as being a simple consequence of the fact that 
the Fermi surface does not intersect the line of d-wave 
gap nodes, so there are no gapless quasiparticles. 

To obtain a more complete phase diagram, we have 
used the insights obtained from the present exact diag- 
onalization studies. It is, of course, not clear to what 
extent the results from small system studies can be ex- 
trapolated to the N — > oo limit. However, at least in 
the case of the pair-binding energy, since the results ap- 
pear to be relatively insensitive to boundary conditions 
for t' < 0.8t, we feel that we can use these results as the 
basis of a set of plausible conjectures. These arc shown 
as the dashed lines in Fig. [H] Notice that we have su- 
perposed the conjectured phase diagram on the contour 
plot of the pair-binding energy from Fig. [2Ja). Where 
the pair-binding energy on the 4x4 system is large and 
insensitive to boundary conditions, we feel that we are 
on sound grounds when we speculate that this corre- 
sponds to a well developed, nodcless d-wave supercon- 
ducting state in the N — > oo limit. 

Unfortunately, where A p b is small and/or sensitive to 
boundary conditions, this could mean that in the A*" — > oo 
limit, the system has entered a gapless phase, i.e. a FL or 
a nodal d-wave superconductor. However, it could sim- 
ply mean that the minimum gap of a nodeless d-wave su- 
perconductor is small and the corresponding correlation 
length is long. In drawing our conjectural phase dia- 
gram, we have assumed the former, and so, to the extent 
possible, we have drawn phase boundaries between the 
nodeless d-wave phase and various gapless phases along 



contours separating the region of "large" and boundary 
condition insensitive pair-binding energy, to regions with 
smaller, and/or strongly boundary condition dependent 
pair-binding. Clearly, the upper portion of the phase 
diagram [if > 0.8£) is the most speculative portion, in- 
cluding the entire region in which nodal d-wave super- 
conductivity occurs. 

The existence of a tetracritical point in the conjectured 
phase diagram, with the additional implication that there 
exists a nodal d-wave superconducting state for large 
enough t' , follows from the nature of the known phases 
along the edges of the phase diagram. In particular, if 
there is a direct, continuous transition from FL1 to a d- 
wave superconductor, the d-wave superconductor must 
be nodal. However, we cannot rule out the possibility of 
a direct first order transition from FL1 to a d-wave su- 
perconductor, in which case the tetracritical point could 
be replaced by a bicritical point, and the d-wave super- 
conductor could always be nodeless. The portion of the 
phase diagram with t' > 0.8t could also exhibit additional 
broken symmetry phase s 11 ' 12 ' 13 . 



B. The Striped Model 

It is impossible, with any degree of confidence, to use 
the present results to infer anything new about the phase 
diagram of the striped Hubbard model in the N — ► oo 
limit. For t' -C t, the model becomes an array of weakly 
connected two-leg ladders, a problem which was studied 
previously by E. Arrigoni et al. in Ref. [Til . From that 
work, we know that for x < x c ps 0.1, there exists a 
nodeless "d-wave like" superconducting state over a very 
broad range of U/t. In contrast, the oscillatory behavior 
of the pair-binding energy seen in the Fig. [5{c) as t' — > 
is a special feature of the a 4 x 2 ladder, which does not 
occur in a system of weakly coupled, infinitely-long two- 
leg ladders. Combining this observation with the strong 
boundary condition dependence of the pair-binding en- 
ergy apparent in Fig. [21b) (even when the twist in the 
boundary conditions is applied in the direction perpen- 
dicular to the stripe direction), we are forced to conclude 
that the results on the 4x4 stripe lattice are not repre- 
sentative of the N — > oo limit. 

We therefore do not venture to draw even a conjec- 
tural phase diagram for this system. Nevertheless, on 
the basis of the fact that, in Fig. [He), there is an ex- 
tended region with relatively large pair-binding energy 
when t' is a substantial fraction of t, makes plausible the 
speculation made in Ref. [ll| that the nodeless supercon- 
ducting state grows in strength for a substantial range of 
non-infinitessimal t' /t. Furthermore, if indeed there is a 
nodal d-wave state for t' /t near 1, there must also exist a 
Lifshitz-type phase transition to a superconducting state 
with gapless quasi-particlcs at a U/t dependent critical 
value of the inter-ladder couplingii. 
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IV. DISCUSSION 

One issue that is often ignored in discussions of HTC 
is the role of the longer-range Coulomb interactions, d- 
wave pairing avoids the obvious deliterious effects of 
the on-site Hubbard repulsion between electrons, but 
is known^ to be fairly sensitive to longer range repul- 
sive interactions. To address this issue, we have com- 
puted the effect on the pair-binding energy of a nearest- 
neighbor repulsion, V, between electrons. Indeed, we 
always find that the pair-binding energy decreases, more 
or less linearly, with increasing V. However, where the 
pair-binding is strong for V = 0, it remains positive 
up to rather large values of V. For instance, for the 
checkerboard lattice under optimal conditions, U = 8t 
and t' = t/2, A p b is an essentially linear function of V 
which vanishes at V w 1.3t. 

In conclusion, we have studied inhomogeneous Hub- 
bard models, primarily with checkerboard and striped 
patterns, on a 4 x 4 square lattice with periodic bound- 
ary conditions by exact diagonalization. Although the 
existence of the HTC in the uniform Hubbard model is 
still a controversial issued, we have produced clear evi- 
dence that, without considering other non-electronic de- 
grees of freedom such as phonons, strong pairing of elec- 
trons can be achieved from purely repulsive interactions if 



certain modulations of the electronic structure are intro- 
duced. Non-monotonic dependence of the pair-binding 
energy and the pair-field pair-field correlators on the de- 
gree of inhomogeneity (f /t) were found to be generic. 
This observation supports the notion that there is an op- 
timal inhomogeneity for high temperature superconduc- 
tivit y 11 ' 17 ! 18 ' 19 ! 20 ' 21 . Since exact diagonalizaiton studies 
cannot access significantly larger systems, it has not been 
possible to carry out finite size scaling to corroborate this 
conclusion. However, we hope that the present results 
will stimulate further work on larger system using more 
efficient numerical tools such as Density Matrix Renor- 
malization Group (DMRG)^ 2 - or Quantum Monte Carlo 
(QMC). 
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